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ABSTRACT 

The magnetar IE 1547—5408 recently exhibited two periods of outburst, beginning on 2008 October 
3 and 2009 January 22. Here we present an analysis of the persistent radiative evolution and a 
statistical study of the burst properties during the 2009 outburst using the Swift X-ray Telescope 
(XRT). We find that the 1-10 keV flux increased by a factor of ~ 500 and hardened significantly, 
peaking ~ 6 hours after the onset of the outburst. The observed pulsed fraction exhibited an anti- 
correlation with phase-averaged flux. Properties of the several hundred X-ray bursts during the 2009 
outburst were determined and compared to those from other magnetar outburst events. We find that 
the peaks of the bursts occur randomly in phase but that the folded counts that compose the bursts 
exhibit a pulse which is misaligned with the persistent pulse phase. We also report a correlation 
between burst hardness and flux. We compare the hardness-flux evolution of the persistent emission 
of both outbursts to those from other magnetars and find that although there does exist an overall 
trend, the degree of hardening for a given increase in flux is not uniform from source to source. These 
results are discussed in the context of previous results and within the magnetar model. 

Subject headings: stars: neutron — pulsars: individual (IE 1547-5408, PSR J1550-5418, SGR 
J1550— 5418) — X-rays: bursts — X-rays: general 



1. INTRODUCTION 

Originally classified as two distinct types of ob- 
jects, anomalous X-ray pulsars (AXPs) and soft gamma 
repeaters (SGRs) are now generally accepted to be 
highly magnetized neutron stars, magnetars, with mag- 
netic fields B > 10 14 — 10 15 G (for re views, see 
iWoods & Thompson! [200l iMereghettil [2008). Magne- 
tars differ from 'normal' rotation-powered pulsars in 
that their X-ray luminosities exceed the rate of en- 
ergy released from their spin-down. They are believed 
to be powered by the decay of their magnetic fields. 
Another common property of magnetars is that they 
exhibit episodes of violent bursting activity. At one 
point thought to be only observed in SGRs, X-ray 
bursts in AXP s were first detected from IE 1048.1-5937 
in late 2001 (|Gavriil et al.l I2002D . Bursts have since 
been detected in se v eral other AXPs (iKaspi et alj|2003t 
Woods et al.l 120051 iKrimm et all 120061 : iGavriil et al.l 



2009). 

IE 1547—5408 was first d iscovered as an X-ra y 
source by the Einstein satellite (|Lamb fc Markcrt 1981). 
It was identified o nly r ecently as a magnetar by 
IGelfand k, Gaenslen (|2007l ) based on its X-ray spectrum 
and infrared flux, as well as its possible association with 
the supernova remnant G327.24— 0.13 . Pul sed radio 
emission was detected by iCamilo et all (|2007h at a pe- 
riod of ^2 s, the shortest period of all known magne- 
tarfl An XMM-Newton observation in 2007 showed a 
significant flux enhancement fr om what was previously 
measured (lHalpern et al.l I2008T ) . thus revealing that an 
X-ray outburst event had occured between 2006 and 

1 Department of Physics, Rutherford Physics Building, McGill 
University, 3600 University Street, Montreal, Quebec, H3A 2T8, 
Canada 

2 see the McGill SGR/AXP Online Cat alog, 
http: / /www. physics, mcgill.ca/ ~pulsar/magnetar/main.html 



2007. Another outburst event oc cured on 2008 Oc tober 
3 and was detected by the Swift (jlsrael et al.ll2010j ) and 
Fermi (|Kaneko et al.l(2010D satellites. IE 1547-5408 en- 
tered yet another active phase on 2009 January 22 when 
hundreds of bursts were d e tected by Swift, INTEGRAL, 
and Fermi. iTiengo et al.l (|2010f) report the appearance 
of dust scattering rings around IE 1547—5408 after the 
2009 burst in follow-up observations with XMM-Newton 
and Swift. The energy released by the event causing the 
rings was estimated to be 10 44 — 10 45 ergs. From these 
rings, a distanc e to th e source was determined to be 3.9 
kpc. INg et al.1 (l20ll report on Chandra observations 
beginning ~ 2 days following the outburst. They note a 
lack of spectral variation during the flux decay as well as 
an anti-correlation between pulsed fraction and phase- 
averaged flux. 

The origin and flux evolution of magn etar outbursts 
are n ot presently well understood (but see lPerna fc Pons! 
1201 If ), and the sample of observ ed outbursts is still re l- 
atively small (for a review, see iRea fc Espositol l20ll . 
Hence, it is important for each magnetar outburst to be 
promptly observed and studied. The magnetar model 
suggests that outbursts are a result of magnctosphcric 
twists of the magnetic field structur e following some 
form of internal energy o r stress release (| Thompson et al.l 
l2002tlBeloborodovl2009D . which predicts a correlation be- 
tween hardness and X-ray flux in magnetar outbursts. 
Indeed a universal relationship between flux increase 
and spectral hardening might be expected, if not from 
source to source, at least for individual sources. Also, 
the origin and physics of X-ray bursts from magnetars 
are poorly understood. Bur sts have been p roposed to 
be magnetospheric in origin (jLvutikovl [20 03) as well as 
originating from stresses in t he crust of the neutron star 
(Th ompson fc Duncar3ll995h . Detailed statistical stud- 
ies of magnetar bursts can help determine their proper- 
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ties for comparison with model predic tions, yet have only 
been done for three sources thus far (jGogiis et al J 120011 
IGavrhl et alll200l . 

This paper presents the results of an analysis of the 
persistent emission from the 2009 event as well as a sta- 
tistical study of the bursts using observations from Swift. 
The burst study will focus primarily on bursts from the 
2009 outburst since the number (~ 400) is much higher 
than in t he 2008 event, in which 8 bursts were detected 
by Swift ([Israel et alJl2010f ). A summary of the observa- 
tions is presented in §2. The analysis performed on the 
data as well as the results for the persistent emission are 
reported in §3.1. The analysis and results for the burst 
study are presented in §3.2. The results and their possi- 
ble physical interpretations are discussed in §4. Finally, 
our findings are summarized in §5. 

2. OBSERVATIONS 

The data presented in this pap er were obtained us- 
ing the X-Ray Telescope (XRT) (|Burrows et all I2005T ) 
on the Swift satellite. The XRT uses Wolter-I optics 
and an XMM-Newton/EPIC MOS CCD detector to pro- 
vide rapid imaging and spectra of X-ray transients in the 
0.5-10 keV energy range. During the 2008 outburst of 
IE 1547-5408, the Swift Burst Alert Telescope (BAT) 
triggered at 09:28:08 UT on 2008 October 3 and Swift 
promptly slewed to the source position. The XRT began 
taking observations 99 s after the trigger. On 2009 Jan- 
uary 22, the BAT triggered at 01:32:41 UT and the first 
XRT observation began ~50 min later. Table [T] shows a 
summary of the observations following the 2009 outburst 
event as well as two observations preceding the outburst. 
The observations f ollowing the 2 008 o utburst were previ- 
ously presented by llsrael et al.l ()2010[ ): here we focus on 
the 2009 event. 

Cleaned data products in both windowed-timing (WT) 
and photon-counting (PC) modes were obtained from 
the HEASARC Swift Archive. Data were reduced to 
the barycentre using the position of IE 15 47—5408, 
15 h 50 m 54Tl, -54°18'23'. / 7 (jCamilo et al.l2007h . For WT 
mode, a source region consisting of a 40-pixel-long strip 
centered on the pulsar position was extracted. A back- 
ground region of the same size was extracted on a source- 
free position away from IE 1547—5408. For the PC mode 
data, an annulus with outer radius 20 pixels and 4 pixel 
inner radius was extracted for the source, and an annular 
region of 20 pixel inner radius and 50 pixel outer radius 
was used for the background. The 4 pixel inner region 
was excluded to avoid pileup. 

3. ANALYSIS AND RESULTS 

3.1. Persistent flux evolution 

In order to investigate the behavior of the persistent 
flux of IE 1547-5408 during the 2009 outburst, the ob- 
servations were fitted with spectral models. Spectra were 
extracted from the event lists using xselect. Spectral 
fitting was performed using the XSPEC0 package ver- 
sion 12.6. The spectra were grouped with a minimum of 
20 counts per energy bin. Ancillary response files were 
created using the FTOOL xrtmkarf and the standard 
spectral redistribution matrices from the Swift CALDB 

3 http://xspec.gfsc.nasa.gov 



were used. Bursts were then removed from the observa- 
tions using the method described in i)3.2l 

The first two observations following the BAT trig- 
ger were split into shorter intervals ^2 ks in length, 
since the spectral properties were evolving rapidly dur- 
ing that time. These split observations occured during 
the fi rst day of the outb urst when the dust scattering 
rings (|Tiengo et al.H2010T) were not fully resolved by the 
Swift XRT. For these observations, the radiative prop- 
erties of the source cannot be simply disentangled from 
the delayed emission from the dust scattering rings (A. 
Tiengo & P. Esposito, private communication). Appro- 
priate modelling of this is currently under investigation 
(Tiengo et al., in preparation) but beyond the scope of 
our work. Following the first day after the outburst, the 
dust scattering rings were outside of the extraction re- 
gion and so those observations are not affected by dust 
scattering. 

The spectra were fitted with a photoelectrically ab- 
sorbed blackbody with an added power-law component. 
To determine Nh, we fit observations having an exposure 
time greater than 3 ks jointly with a single Nh. The pa- 
rameters kT, r and their normalizations were allowed to 
vary in these fits. This resulted in N H = 3.24(5) x 10 22 
cm~ 2 which is consistent with the value of 3.llg'g x 10 22 
cm~ 2 reported in lGelfand fc Gaenslerl ()2007l ) . though is 
somewhat lo wer than the v alue of 4.1(1) x 10 22 cm" 2 
measured bv lNg etldl (|2011[ ). All the Swift observations 
were subsequently fit with the column density fixed to 
our best-fit v alue. The 2008 da ta have been previously 
presented bv llsrael et a l. (2010), whose results are gen- 
erally consistent with those of our analysis for the same 
time period, so they will not be presented here. Figure 
Q] shows the results of fitting of the observations follow- 
ing the 2009 outburst. The split observations, affected 
by dust scattering, are the first 10 data points in Fig- 
ure [T] following the trigger, and so should be regarded 
with caution. The two observations preceding the out- 
burst were fit with only a blackbody and no power-law 
component. This is because the power-law index could 
not be constrained due to a paucity of counts. The fits 
in general were excellent; the goodness-of-fit statistic % 2 
ranged between 0.61 and 1.48 with a mean of 1.06. 

The peak of the persistent emission occured ~6 hr af- 
ter the BAT first triggered on IE 1547-5408 in 2009 
January. The peak was almost 3 orders of magnitude 
higher in flux than the pre-outburst emission. This was 
accompanied by a hardening of the spectrum in the 1-10 
keV band, as can be seen in the falling spectral index 
and rising kT. After the peak, the spectral index soft- 
ened and kT fell as the flux dropped. To characterize the 
flux decay of the outburst, a power law decay was fit to 
the data following the first day of the outburst. The first 
day was not included as the source decay is superimposed 
by delayed emission from the dust scattering rings. The 
power law decay is descibed by F = A(t — to) a , where F 
is the unabsorbed flux, A is the normalization, a is the 
power-law index, and to is the time of the BAT trigger. 
The decay was described by a power-law with an index 
of -0.24 ±0.02. Although the xll" of the fit was 4.1/23, 
a power-law model fit the data much better than an ex- 
ponential decay. To determine the blackbody radius of 
the emitting region as shown in Figure [TJ a distance of 
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3.91 ± 0.07 kpc from lTiengo"etaLl (f2010h was assumed. 

To measure pulsed fractions and fluxes, the burst- 
removed time series (see §3.2) were folded at the rota- 
tional ephemeris derived from contemporaneo us RXTE 
observations of IE 1547-5408 (|Dib et al.ll20TTI ). Specifi- 
cally, we used spin frequency v — 0.48259615(3) Hz, with 
frequency derivative v = —5.12(2) x 10~ 12 s~ 2 at refer- 
ence epoch MJD 54854.0 or 2009 January 23. The pulsed 
flux was calcu lated using an R MS method according to 
the formula in lDib et alT(|2008D . with 7 har monies. The 
pulsed fraction was determined by dividing the pulsed 
flux by the phase-averaged flux. A pulsed flux and frac- 
tion were measured for each WT observation with an 
exposure time > 3 ks for which the ephemeris was valid. 
As seen in Figure [TJ the pulsed fraction decreased as the 
flux increased. After ~40 days, the pulsed fraction had 
not yet recovered to its pre-burst level of ^0.25. 

The pulsed flux evolution around both the 2008 and 
2009 outburst events is presented in Figure [5J The 
Swift data confirm the evolution that is ob served in 
the R XTE data whic h are presented b riefly in INg et all 
(|2011[) and in detail in lDib et all (|2011l ). The pulsed flux 
enhancement ^11 days following t he initial trigger of the 
2008 event, although not noted bv llsrael et al l (|2010t ). is 
clearly present in the Swift data. Thus, the 2009 event 
showed a much smaller increase in pulsed flux than did 
the 2008 event, while the opposite is true of the total 
flux. 

Figure [3] shows an anti-correlation between pulsed frac- 
tion and unabsorbed flux. It includes o bservatio ns from 
both the 2008 and 2009 outbursts. INg et all (|2011[ ) 
present this trend using Chandra and XMM data as well 
as the burst-removed Swift data presented in this paper. 

3.2. X-Ray Bursts 

Bursts w ere identified i n a m anner similar to that de- 
scribed in iGavriil et al.l (|2004l ). The event lists were 
binned at a 1/16 s time resolution and a mean number 
of counts per bin was calculated for each Good Timing 
Interval (GTI). The number of counts in each bin, rii, 
was compared to the GTI mean counts, A, according to 
the probability Pi of rii occuring randomly, 



Time bins that had P t < 0.01/N, where N is the total 
number of time bins in the GTI being searched, were 
identified as part of a burst. Since the mean number of 
counts in a GTI can be overestimated due to contam- 
ination from bursts in other time bins, the procedure 
above was repeated iteratively, each time removing the 
bins that were identified as containing bursts, until no 
further bursts were identified. The above procedure was 
then repeated for 1/32 s and 1/64 s time resolutions to 
improve sensitivity to bursts of different durations. 

Since several bins identified with bursts can be part of 
the same burst, a burst was defined by its peak. The peak 
of a burst was determined by first finding the minimum 
time to accumulate 10 counts, using unbinned event data. 
The midpoint of the time spanned by these 10 counts was 
defined as the peak. A search for a peak was done within 
0.5 s on each side of an identified bin. This definition of 
burst peak is independent of binning and, for bursts with 



durations shorter than 1 s, will merge all of the identified 
bins into a single burst. Bursts within 1 s of each other 
that were merged into a single burst were identified when 
selecting the background (see below) and their properties 
were measured separately. 

Once identified, to remove the bursts from the event 
lists, the lists were divided into full periods of the pulsar, 
starting with the first event of the observation as a refer- 
ence point. The periods that contained bursts were iden- 
tified and all counts which arrived in that interval were 
removed from the event list. This was done to ensure 
equal exposure to all pulse phases in the pulse profile. 

3.2.1. Burst Statistics 

For each burst, a fluence, Tgo, rise time (t r ), and fall 
time (tf) we r e me asured with the same analysis as in 
IGavriil et all((200l . using the unbinned event data. The 
fluence was determined by first measuring a background 
count rate in hand-picked regions on either side of the 
burst. The background region by default was between 1 s 
and 2 s from the burst peak to either side of the burst, but 
was manually adjusted for most bursts in order to avoid 
contamination from other nearby bursts. The cumulative 
background-subtracted counts were then fit with a step 
function, using data point from the hand-picked back- 
ground region. The height of the step function in counts 
corresponds to the fluence of the burst. The dura- 
tion of the burst is the time between when 5% and 95% 
of the fluence has been accumulated. The burst rise and 
fall times were determined using a maximum likelihood 
fit to a piecewise function with an exponential rise and 
an exponential decay. Peak fluxes in counts per second 
were determined by passing a 62.5 ms boxcar integrator 
through a 250 ms interval (4 boxcar widths) centered on 
the burst peak. The highest count rate measured by the 
integrator was defined to be the peak flux. 

In total, for the 2009 outburst, 424 bursts were iden- 
tified in 86 ks of observations from 2009 January 22 to 
2009 September 30 using 1/16 s time resolution. Thir- 
teen additional bursts were identfied using the 1/32 s 
and 1/64 s time resolutions for a total of 437 bursts. Of 
those identified, 34 had too few counts to permit a reli- 
able measure of fluence, and for 32, the exponential rise 
and decay were not successfully fit. For 64 of the bursts, 
properties could not be measured because they were too 
close to another burst to allow a reliable background es- 
timate between them. Four bursts were too close to the 
edge of a GTI to calculate a background count rate. In 
summary, 303 bursts could be fully analysed and only 
these 2009 bursts will be considered henceforth. Exam- 
ples of different bursts are shown in Figure |4j Only two 
bursts were found in the XRT observations of the 2008 
event; their properties were similar to those during the 
2009 outburst. 

Figure [5] shows the distributions of the burst proper- 
ties, grouped in logarithmic bins. The Tgo, rise time, fall 
time, and the ratio of rise to fall time distributions have 
been fit with log-normal distributions using maximum 
likelihood fitting. The Tgo distribution has a mean of 
305 ms and a range for one standard deviation of 140 - 
662 ms. For the rise time distribution, we find a mean of 
39 ms and a range for one standard deviation of 14 - 109 
ms. For the distribution of fall times, we find a mean 
of 66 ms and a range of 24 - 182 ms for one standard 
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deviation. The mean of the t r /tt distribution is 0.59 
with a range for one standard deviation of 0.21 - 1.66. A 
summary of the measured quantities is provided in Ta- 
ble 2. There we also provide these quantities, when avail- 
able, for the three other sources for which such statistical 
analyses have been done. We note in particular that the 
average burst duration for IE 1547—5408 is longer than 
for all others measured thus far. This result is further 
discussed in §4.2. 

Figure |5] shows the fiuence and peak flux distributions, 
in counts and counts per second, respectively, which 
have been fitted with a power law using a least-squares 
method. For low fluences or peak fluxes, the number 
is underestimated since the burst detection algorithm is 
less sensitive to these bursts. Therefore, some of the 
points with low fiuence or peak flux were not included 
in the fits. The best-fit power-law index for the fiuence 
distribution is —0.6 ±0.1. For the peak flux distribution, 
the power-law index is —0.34 ± 0.12. This index is quite 
shallow and, despite the omission of the first two points 
in the fit, may still be affected by the bias in the burst 
search. 

For other AXPs, t here is evidence that burs ts tend to 
arrive on pulse (e.g. IGavriil et all 120021 12004ft . For IE 
1547—5408 this does not seem to be the case, as burst 
peaks are distributed randomly in phase (see Fig. [Th). 
This is simila r to what is observ ed in SGRs 1806—20 
and 1900+14 (jPalmerl [l999l [20021 ) . However, when the 
individual photon arrival times that are part of bursts 
are folded, a strong pulse is observed. The peak of this 
'pulse', shown in Figure [7)d, is not aligned with the peak 
of the quiescent pulse profile. Figure [7^ presents this 
quiescent pulse profile obtained using the burst-removed 
data. Since the pulsed fraction of the first two obser- 
vations following the 2009 BAT trigger is significantly 
lower and the persistent flux is significantly higher than 
in the subsequent observations, including them in the 
profile reduces the pulse amplitude. Thus, in order to 
better show the persistent pulse profile, Figure [7^ does 
not include those first two observations. In order to de- 
termine whether the burst count pulse was dominated by 
a handful of bright bursts, the 15 brightest bursts were 
removed and the counts were refolded. The profile did 
not change significantly and still displayed a pulse at a 
similar phase. 

Panels c and d of Figure [7] separate the symmetric 
bursts from slow-fall bursts. Symmetric bursts were de- 
fined as bursts with 0.5 < t r /tf < 2. Bursts with tf 
greater than 2 t r were defined as slow-fall bursts. Of the 
303 well measured bursts, 158 were classified as symmet- 
ric and 116 as slow falls. The remainder, 29 bursts, were 
those with rise times greater than twice their fall times. 
Curiously, the folded slow-fall burst counts exhibit a 
much stronger pulse than do the folded symmetric burst 
counts, with a xt — 43 for the null hypothesis, com- 
pared to xt = 6 f° r the folded symmetric bursts. This 
symmetric/slow- fall definition is somewhat arbitrary and 
these two classes of bursts are by no means distinct pop- 
ulations. We use this distinction only to demonstrate 
that the more symmetric bursts tend to be more pulsed 
than those that are less symmetric. 

3.2.2. Burst Spectroscopy 



Spectra within the Tgo interval of each burst peak were 
extracted and grouped with 20 counts per bin. Back- 
ground spectra were taken from 1 min on either side of 
the burst peak extracted from the burst-removed event 
data. The 46 bursts with fluences over 200 background- 
subtracted counts were fitted with a photoelectrically ab- 
sorbed power law with Nh fixed to 3.24 x 10 22 cm" 2 
as measured from the fit to the persistent emission (see 
§3.1). The spectra were fitted using XSPEC. The mean 
of the measured spectral indices of the bursts was found 
to be r = 0.17, with a standard deviation of 0.33, where 
N(E) cx E~ r . Fitting with more complicated models 
was attempted, but parameters could not be successfully 
constrained because of the low number of counts. As 
shown in Table 2, the average T we measure is signifi- 
cantly harder than that measured in the only other mag- 
netar outburst for which the value is reported. This is 
discussed further in §4.2. 

In order to determine the effect of pileup on the bursts, 
power-law spectra were fit to the two bursts with the 
highest peak flux while excluding a central region ranging 
from 1 to 15 pixels in radius. The power-law indices from 
the fits were plotted ag ainst radius of exclu sion region as 
per the prescription of lRomano et al.l (|2006l ). The power- 
law index did not change significantly as a function of 
exclusion radius. Hence we concluded that the bursts 
were not significantly affected by pileup. 

In order to probe the relation between hardness and 
fiuence for the bursts, hardness ratios were determined 
for each burst by measuring fluences in the 0.5-4 keV 
and 4-10 keV bands. No significant correlation between 
0.5-4 keV/4-10 keV hardness and fiuence was observed 
for bursts from IE 1547—5408. For the 46 most fluent 
bursts, there was also no observed correlation between 
the 0.5-4 keV/4-10 keV hardness or the spectral index 
of the power-law fit and the fiuence. The lack of a detec- 
tion of the correlation between hardness and fiuence in 
the Swift data may be due to the limited spectral range 
of the XRT (0.5-10 keV). However, an anti-correlation 
between power-law index and average flux over Tgo from 
the spectral fits was observed (Fig. [S|). This is consistent 
with a correlation between the hardness and the magni- 
tude of a burst. 

Spectral features in magnet ar bursts have been 
repor te d for some so urces (iStrohmaver fc Ibrahim! 
20051 llbrahim et all 12001 IGavriil et al.l 120021: 
Ibrahim et al.l 120031: iWoods et all 120051 : IGavriil et all 
20091 : IKumar fc Safi-Harbl 12010ft . although most of 
the features have been discovered above 10 keV. The 
individual burst spectra for IE 1547—5408 showed 
no evidence for any spectral features in the 1-10 keV 
range. Since spectral features may be masked by the 
low count rates in the individual spectra, we combined 
the burst spectra into a single average spectrum. The 
burst and background spectra and response filestend to 
bere combined using the FTOOL addspec. Spectral 
bins were grouped with a minimum of 20 counts per 
bin. The average spectrum was fit with a blackbody, 
a power law, a blackbody with an added power law, 
and a Comptonized blackbody; all four models were 
photoelectrically absorbed. All of the models gave 
acceptable fits with xt ~ lj thus, no significant spectral 
features were observed in any of the residuals. 

Fluences and peak fluxes in counts and counts per sec- 
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ond, respectively, were converted into cgs units using the 
results of the spectral fits. The fluxes from the power- 
law fits were multiplied by the burst durations and com- 
pared to the fluences measured in counts. The propor- 
tionality constant between the two was determined to be 
2.23 x 10 -10 ergs cm" 2 counts -1 . In using this single 
conversions factor between counts and ergs cm -2 we are 
effectively assuming an average spectral model for the 
bursts which is not generally true. The conversion from 
counts to cgs units in reality is different from burst to 
burst depending on their spectra. Values in cgs units, 
using this single conversion factor, are indicated on the 
top axes of the fluence and peak flux distributions in 
Figure El 

4. DISCUSSION 

We have presented Swift XRT observations of the 2008 
and 2009 outburst events of IE 1547-5408. In partic- 
ular, we considered the behavior of this source's persis- 
tent flux in its 2009 outburst, which was qualitatively 
different from that in 2008. Specifically, the 2009 event 
showed a much smaller increase in pulsed flux than in 
2008, while the opposite is true of the total flux. In this 
work, we have shown that for the 2009 outburst, the 
source spectrum became harder (though the exact de- 
gree of hardening is presently difficult to know, due to 
the contaminating dust scattering) as the flux increased 
after the initial trigger, and we found an anti-correlation 
between the pulsed fraction and the flux. 

We have also presented a detailed study of the sev- 
eral hundred X-ray bursts detected by Swift following 
the 2009 event. Noteworthy results we have found in- 
clude a correlation between burst flux and hardness. We 
also note that burst counts tend to contribute very sig- 
nificantly to the pulsed fraction, even if the phases of 
burst peaks are randomly distributed in pulse phase. 

Next we discuss our results in the context of the mag- 
netar model, as well as in comparison with those for other 
similar outbursts. 

4.1. Persistent Flux 

Figure [T] shows that as the X-ray flux increased, the 
blackbody temperature increased and the power-law in- 
dex decreased during the 2009 outburst. Such a hard- 
ness/flux correlatio n was seen also in the 2008 outburst 
(jlsrael et al.l 120101 ) and ubiquitously in other magne- 
tar outbursts (e g. IKaspi et al.ll2003t IWoods et afl 120041: 
llsrael et al.ll2007D . In the case of IE 1547-5408, however, 
we regard the precise values of kT and V as shown in Fig- 
ure [T] on the first day after the outburst with caution, as 
they are very likely biased by dust scattering. The hard- 
ness/flux correlation for this outburst on the day of the 
trigger is likely to be robust, since dust scattering softens 
photons. However, note that if the dust scattering rings 
are caused by a burst that is much harder than, and sep- 
arate from, the persistent emission, the dust scattering 
could effectively cause the overall emission to harden in 
relation to the true persistent emission. 

The thermal emissio n of magnetars, in the twisted 
magnetosphere model ()Thompson et al. 2002), has as 
its origin heating from within the star, due to the de- 
cay of the strong internal magnetic field. The resulting 
thermal surface photons are thought to be scattered by 
currents in the atmosphere, resulting in a Comptonized 



blackbody- like spectrum whi ch is often modelled with a 
blackbody plus a power law (jLvutikov fe Gavriii| [2006'). 
The magnetospheric currents are prese nt due to 'twists' 
in th e field structure, either global (I Thompson et ah 
2002), or, more likely, in localized regions (jBeloborodov 
2009). In either case, the current strength, hence de- 
gree of scattering, increases with increasing twist magni- 
tude. Moreover, return currents provide an additional, 
external source of surface heating in addition to the in- 
ternal source. The increase in flux during a magnetar 
outburst is thus theorized to be caused by an internal 
heat-releasing event that may significantly increase the 
surfa ce temperature (see, for example, IQzel fc Guverl 
2007), further twist the magnetospheric field, and in- 
crease the external return-current heating. Thus, a corre- 
lation between hardness and flux is generically expected 
in the twisted- magnetosphere model (jLvutikov fc Gavriill 
2006), in agreement with observations of IE 1547—5408 
and other magnetars. 

It is interesting to compare the hardness/flux corre- 
lations seen for IE 1547—5408 with those observed for 
other magnetars to see if there exists a universal rela- 
tionship between increase in X-ray flux and hardness, as 
might be expected in the above scenario. Figure [9] shows 
the fractional increase in the 4-10/2-4 keV hardness ra- 
tio determined from fluxes as a function of 2-10 keV 
unabsorbed flux fractional increase for six different mag- 
netar outbursts in four different magnetars. We chose 
to consider a flux hardness ratio as a measure of hard- 
ness (rather than e.g. power-law index) as it is defined 
independent of spectral model, and also is instrument- 
independent. In order to determine the hardness of each 
source, spectral model parameters from the literature 
were compiled and input to XSPEC. The 4-10 keV and 
2-4 keV fluxes were then determined using the XSPEC 
flux command and the hardness was calculated from 
their ratio. For references that did not report a 2-10 
keV unabsorbed flux, one was determined using XSPEC 
using the reported model paramenters. As is clear from 
the Figure, an overall trend is apparent, although there 
is significant variation from source to source. For ex- 
ample, for similar flux enhancements over the quies- 
cent level, IE 2259+586 became much harder than did 
IE 1048.1-5937, with the behavior of IE 1547-5408 
lying somewhere in between. Thus, though Figure [5] 
demonstrates that a hardness/flux correlation is indeed 
generically observed, there does not appear to exist a 
universal law linking the degree of flux increase over the 
quiescent level with the degree of flux hardening. This 
observation will require accommodation in any detailed 
model of magnetar outburst emission. 

As shown in Figure G2 the flux of IE 1547-5408 was 
anti-correlated with the pulsed fraction. While immedi- 
ately preceding the 2009 outburst the pulsed fraction was 
over 20% (Fig. [IJ, immediately afterward, it dropped 
to nearly negligible levels in the 1-10 keV range, likely 
due to the emission being dominated by the scattered 
emission from the dust scattering rings. However, fol- 
lowing the first day, the pulsed fraction was still signif- 
icantly lower than the preburst pulsed fraction. This 
anti-correlation between pulsed fraction and flux is pre- 
sented with the additio n of Chandra and XMM-Newton 
data in INg et al.l (|2011l ). who suggest it could be due to 
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an increase in the emitting area of a thermal hot spot, 
which would cause the emission to be observable during a 
greater portion of the phase of the pulsar. We note that 
both correlations (iGotthelf et al.l|2004t [Zhu et al.1 1 20081) 
and anti-correlations ([Israel et al.ll2007t iTam et aljboOSl i 
between pulsed fraction and flux have been observed in 
magnetar outbursts. Different behaviors can result de- 
pending on the location of the emission region and the 
viewing geometry. Although detailed modeling of spec- 
tral, pulse morphology and pulsed flux changes during 
magnetar outbursts is beyond the scope of this paper, 
recent attempts at unified modelling of magnetar sur- 
face and magneto sphere emission geo metries and emis- 
sion mechanisms (Alb ano et al.|[2010l ) could in principle 
use data like ours to constrain the hot spot location and 
size. This seems to us to be a good future avenue for 
investigation. 

Interestingly, iKaneko et al.l (|2010f ) report a —55% 
pulsed fraction at —100 keV as observed by Fermi GBM 
—30-40 minutes before the BAT trigger. They report 
that this is the first time that pulsations unrelated to a 
giant flare have been observed at -100 keV for an SGRE|. 
However, AXPs have previously been shown to exhibit 
pulsations at —100 keV with pu lsed fractions as high 
as 100% (e.g. iKuiper et all 12006ft . Thus it is possible 
that the measured 55% actually decreased from a higher 
pulsed fraction leading up to the outburst event, i.e. the 
—100 keV pulsed fraction may have behaved similarly to 
that in the 1-10 keV band. Future hard X-ray telescopes 
with focusing optics, such as NuSTAR, will allow much 
easier measurements of the pulsed fraction of magnetars 
at high X-ray energies. 

4.2. Bursts 

Previous studies of SGR bursts have shown that their 
energies, and thus fluences, follow a power- law distribu- 
tion d N/dE oc E~ a with a equa l to -5/3 ()Cheng et al.l 
[19961: IGoeiis et al.l Il999l |2000). It has been noted 
that this is similar to the Gutenberg-Richter law for 
earthquakes and to energy distrib utions of solar flares 
(|Crosbv et al.lfT99llLu et al.lU99l . Both IE 1547-5408 
and IE 2259+586 also follow this distribution. In this 
work, the fluence distribution of IE 1547—5408 is found 
to have a power-law index of —0.6 ± 0. 1 whic h corre- 
sponds to dN/dE cx E~ 16 . IGavriil et all (|200l find an 
a of 1.7 ± 0.1 for IE 2259+586, similar to the values 
found for SGRs, further reinforcing the similarity in this 
particular behavior among AXPs and SGRs. 

For the 2009 outburst, INTEGRAL observations which 
cover an energy range of > 80 keV, show different burst 
properties fro m those determin e d usin g the Swift XRT 
observations. Sa vchenko et al.1 (|2010l ) report a 68-ms 
mean duration derived from a log-normal distribution 
with a scatter of 30 - 155 ms. This is much shorter than 
the 305-ms duration determined in this work (see Ta- 
ble [2]) . This discrepency may be due to the difference 
in energy coverage; perhaps bursts have different mor- 
phologies at differe nt energies. How e ver, t he definition 
of duration used bv lSavchenko et a"H (120101) d i ffers f rom 
the Tg definition used here. I Savchenko et all (|2010[) de- 
fine the burst duration as the time between the moment 

4 Note that [Kancko ct al. (2010) use the designation SGR 
1550-5418 for IE 1547-5408. 



when the count rate rises above 5a to when the count 
rate drops below 3a. When applying their definition of 
duration to the bursts identified in our study, we find, 
for a log-normal distribution, a mean duration of 101 ms 
and a range for one standard deviation of 59 - 173 ms, 
closer to but still somewhat longer than their measure- 
ment, suggesting a possible energy dependence of burst 
duration. We do not, however, detect any significant 
difference in the durations measured using 0.5-10 keV 
counts with those measured using 2-10 keV counts. 

The properties of bursts from the 2009 outburst event 
of IE 1547—5408 are re miniscent of those f rom the out- 
burst of IE 2259+586 (|Gavriil et al.l 120041) . Both have 
a significant number of short spikes like those found 
in SGRs, and a set of bursts with long pulsating tails 
like those found in burst studies of other AXPs. Al- 
though we do not find any bu rsts with long pu l sating 
tail s in the XRT observati ons, ISavchcn ko et al.l (2010) 
and lMereghetti et al.l ((200S) find two such bursts. These 
were not found in our analysis because Swift was not 
observing IE 1547-5408 when they occured. Table [2] 
compares the properties of bursts from IE 1547—5408 
to those from outbursts from other magnetars. We note 
that the average durations of bursts from IE 1547—5408 
appear to be longer than for those in other sources. How- 
ever, this could be an artifact as the burst properties 
for the other tabulated sources were determined with 
RXTE. The larger collecting area of RXTE allows it to 
detect bursts that are fainter than those Swift can de- 
tect. If faint, short bursts are missed by the XRT, the 
mean burst duration (as well as t r and tf) may be over- 
estimated. Also, the energy range of RXTE (2-60 keV) 
probes higher energies and so differences may be due to 
the energy dependance of burst properties. A detailed 
statistical study of IE 1547-5408 bursts with RXTE is 
needed to clarify this point . 

Although Swift XRT is not ideal for probing the spec- 
tra of magnetar bursts, which have significant flux above 
10 keV, we were still able to draw the following con- 
clusions from our analysis. While we do not observe 
a hardness-fluence correlation for IE 1547—5408, this 
could be due to our limited energy range. Indeed, there 
is a hint of a correlation between T and fluence, how- 
ever is is not statistically significant. Note however that 
we do observe a signifi cant r-flux correlation (Fig. [8]). 
iSavchenko et al.1 (|2010l ). using INTEGRAL data of the 
2009 outburst, find a correlation between burst hardness 
and count rate. Their hardness ratio is defined as the 
ratio between the Anti- Coincidence Shield (ACS) flux, 
which is sensitive to photons above 80 keV, and a 20-60 
keV flux from the ISGRI instrument. This is also con- 
sistent with a correlation bet ween hardness and burst 
magnitude for IE 1547-5408. IGavriil et all (|2004l ) also 
find a correlation between hardness and fluence for IE 
2259+586. For the SGRs, on the other hand, an anti- 
correla tion between hardness a nd fluence has been ob- 
served (jGogiis et al.lll999ll2000l) . Table [2] also shows that 
the IE 1547—5408 bursts from Swift are much harder 
than those from IE 2259+586. In light of the observed 
hardness-fluence correlation in AXP bursts, this is not 
a suprising result. The 28 most fluent bursts from 
IE 2259+586 for which spectral indices w ere measured, 
have fluences of — 10 -9 — 10 -8 erg cm -2 (|Gavriil et al.l 
2004). This is to be compared with the 46 most fluent 
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bursts for which T was measured here for IE 1547—5408, 
that span a fluence range of ~ 10~ 8 — 10~ 7 erg cm -2 . 
This may account for the harder average spectral index 
for bursts from IE 1547-5408. 

In the magnetar model, two mechanisms have 
been suggested for produ cing magnetar bursts. 
iThompson fc Duncan! (|1995j ) suggest that stresses 
due to the strong magnetic fields present inside magne- 
tars are able to crack the crust of the neutron star. This 
cracking releases a plasma fireball into the magneto- 
sphere. The strong magnetic fields can hold the fireball 
above the fracture site. The suspended fireball can heat 
the surface thus causing an extended cooling tail. Since 
the strongest surface fields are located near the polar 
caps, these fracture events would occur preferentially 
near the poles which would result in an observed phase 
de pendence of the b ursts. Another mechanism, proposed 
bv iLvutikovl ((2003), suggests that bursts are caused by 
reconnection events initiated by the development of a 
tearing mode instability in the magnetically dominated 
relativistic plasma in the magnetosphere. In this case, 
bursts occur randomly in phase. This mechanism 
should also produce shor ter more symm etric bursts 
than in crustal fracture. ILvutikovl (|2003[ ) states that 
the hardness-fluence anticorrelation found in SGRs is 
consistent with magnetic reconnection. For the surface- 
co oling model, the opposite is expected. As discussed 
by ILvutikovl (|2003[ ) , both mechanisms could easily be at 
w ork. 

iWoods et al.1 (j2005[ ) suggested that observationally, 
there are two types of magnetar bursts. Type A bursts 
are nearly symmetric and are typical of SGR bursts. 
Type B bursts have been observed in AXPs and are char- 
acterized by a short spike followed by a long tail, typi- 
cally much longer than the rotational period of the pul- 
sar. Pulsations have been observed in these tails. Since 
the Type B bursts observed in AXPs occur preferentially 
in phase and exhibit a hardness-fluence correlation, and 
the Type A bursts in SGRs are distributed randomly 
in phase and have a hardness-fluence anti-correlation, 
IWoods et all (|2005| ) suggest the magnetic reconnection 
mechanism for type A bursts and the surface-cooling 
model for Type B bursts. 

For IE 1547-5408, none of the bursts detected by Swift 
XRT could be classified as Type B. However, the two 
bursts w ith long pulsatin g tails foun d in INTEGRAL 
data by iSavchenko et al.l (|2010l ) and iMereghetti et al.l 
(2009) have typical Type B morphology. For the bursts 
in this work, over half of the bursts were classified as 
symmetric, which nominally corresponds to Type A (see 
top- right panel of Figure [4] for example). About 40% of 
the bursts were classified as slow-fall bursts, which are 
actually closer to Type A in morphology than Type B, 
although they are not very symmetric. Thus, the bursts 
from IE 1547—5408 are n ot easily class i fied in to Types 
A and B as described by IWoods et all (|2005l) . More- 
over, for both the symmetric and slow-fall bursts, the 
folded photon arrival times exhibited a clear phase de- 
pendence, although for the slow-fall bursts this 'pulse' is 
much stronger. That the symmetric bursts show some 
pulse modulation is suggestive of a difference with 'clas- 
sical' Type A bursts, or that pulse phase analyses of the 
latter should be attempted using all burst counts, as they 
too may be pulsed. That the burst counts are pulsed 



indicates that burst emission comes from a preferred re- 
gion in rotational phase, be it on or near the surface or 
high in the magnetosphere, even if the burst peaks arrive 
randomly in phase. The offset between the burst counts 
pulse peak and the persistent pulse peak seen when com- 
paring Figures 7b and 7e demonstrates that the preferred 
burst emission region has location and geometry similar 
to, but distinct from, that producing the persistent pul- 
sations. 

The total energy released in the Swz/t-detected bursts 
was ~ 1 x 10 40 erg in the 1-10 keV range. This is 
much lower than the 1-10 keV energy released from the 
persistent emission between 2009 January 22 and 2009 
September 30 of ~ 9 x 10 41 erg. For reference, the en- 
ergy released fr om the spin-down in that same period is 
~ 5 x 10 40 erg. IWoods et alJ (pOOl find that for SGRs, 
the energy released in the bursts is higher than that re- 
leased in the persistent emission, but for IE 2259+586, 
the opposite is true. In this regard, the 2009 outburst of 
IE 1547-5408 is more like that of AXP IE 2259+586. 

5. CONCLUSIONS 

We have presented an analysis of the persistent 
radiative evolution of the 2009 January outburst of 
IE 1547-5408 from Swift XRT observations. We found 
that ~ 6 hr after the 2009 BAT trigger the observed 
persistent 1-10 keV unabsorbed flux reached a peak of 
~8x 10 -9 ergs cm -2 s _1 , an increase of more than 500 
times the quiescent flux leading up to the outburst. This 
flux evolution is not due solely to the source; in the first 
day, there is also emission from dust scattering rings that 
is delayed emission from an energtic event near the onset 
of the outburst. There was significant spectral harden- 
ing at the outburst as seen in other magnetar outbursts. 
Note tha t the a bsence of spectral variation reported by 
iNg et al.1 (|2011l ) is consistent with our results, as they 
missed the bulk of the spectal changes which occured in 
the first day of the outburst and their Chandra observa- 
tions did not begin until the next day. The pulsed frac- 
tion showed an anti-correlation with the phase-averaged 
flux for both the previous 2008 and 2009 outbursts, with 
both sets of data following the same trend. We have com- 
piled data from this and five other magnetar outbursts 
for four different sources and find a generic X-ray hard- 
ness/flux correlation overall, but with no clear universal 
quantitative relationship between the two. 

We have also presented a detailed statistical analy- 
sis of the several hundred bursts detected during the 
2009 event. The bursts do not easily fall into the 
Type A/Type B classification put forth by IWoods et al.1 
(2005). We found that the peaks of the bursts were ran- 
domly distributed in pulse phase, but that when the in- 
dividual photon counts were folded at the pulsar emphe- 
meris, a clear pulse was present. This phase dependence 
is stronger for those bursts with longer decays than rise 
times than for those bursts that are symmetric. We also 
report a correlation between burst hardness and burst 
flux. 

In many ways, these observations yield more questions 
than answers. The range of observed phenomenology 
in magnetar outbursts seems to increase with each 
event, with few overall trends to assist in constraining 
models emerging. Nevertheless given the paucity of 
events studied in detail, we remain hopeful that through 



perserverence, eventually physical insight will emerge. 
The fast response of telescopes like Swift is crucial to 
this endeavor. For IE 1547—5408, it allowed an analysis 
of the first day of the 2009 event. This is important 
as both the most significant spectral changes and the 
majority of the bursts occured within this period. This 
highlights the necessity of prompt response to magnetar 
outbursts in understanding their nature. 

After submissi on of our manuscr i pt, we became aware 
of the work by iBernardini et all (|2011h on the 2009 
outburst of IE 1547-5408 using Swift XRT data. A 
preliminary comparison of their results with ours shows 
general agreement. In particular their reported overall 
source behavior, namely spectral hardening correlated 
with flux, as well as pulsed fraction anti-correlated with 
flux and blackbody radius, are in broad agreement with 
our results. We further note that they omitted reporting 
on data taken with Swift XRT in the 24 hrs following 
the 2009 trigger, due to the difficulty in handling the 
dust-scattered emission (P. Esposito, A. Tiengo, private 
communication) . 
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Table 1 

Summary of Swift XRT observations of the 2009 outburst of IE 1547-5408 



Sequence Mode Observation date MJD Exposure time Time since trigger 









(TDB) 


(ks) 


(days) 


00090007024 


WT 


2009-01-04 


54835.0 


3.3 


-18.033 


00090007025 


WT 


2009-01-12 


54844.0 


4.2 


-9.066 


00340573000 


WT 


2009-01-22 


54853.1 


6.1 


0.035 


00340573001 


WT 


2009-01-22 


54853.4 


9.2 


0.321 


00340923000 


PC 


2009-01-23 


54854.6 


1.7 


1.575 


00090007026 


WT 


2009-01-23 


54854.7 


8.2 


1.657 


00340986000 


PC 


2009-01-24 


54855.2 


2.9 


2.110 


00030956031 


PC 


2009-01-24 


54855.3 


2.5 


2.195 


00090007027 


WT 


2009-01-25 


54856.1 


3.3 


2.990 


00341055000 


PC 


2009-01-25 


54856.1 


4.0 


3.065 


00341114000 


PC 


2009-01-25 


54856.9 


4.6 


3.851 


00090007028 


WT 


2009-01-26 


54857.1 


3.5 


3.987 


00030956032 


PC 


2009-01-27 


54858.1 


6.2 


4.996 


00090007029 


WT 


2009-01-27 


54858.4 


1.8 


5.330 


00030956033 


PC 


2009-01-28 


54859.2 


5.1 


6.143 


00090007030 


WT 


2009-01-28 


54859.9 


1.9 


6.811 


00030956034 


PC 


2009-01-29 


54860.0 


5.9 


6.941 


00090007031 


WT 


2009-01-29 


54860.7 


2.2 


7.614 


00090007032 


WT 


2009-01-30 


54861.5 


2.9 


8.476 


00030956035 


WT 


2009-01-30 


54861.7 


3.0 


8.677 


00030956036 


WT 


2009-01-31 


54862.1 


3.0 


9.076 


00090007033 


WT 


2009-01-31 


54862.8 


2.5 


9.687 


00030956037 


WT 


2009-02-01 


54863.5 


2.0 


10.476 


00090007035 


WT 


2009-02-02 


54864.6 


1.1 


11.548 


00030956038 


PC 


2009-02-03 


54865.5 


5.9 


12.424 


00030956039 


PC 


2009-02-04 


54866.7 


6.1 


13.637 


00030956040 


WT 


2009-02-05 


54867.2 


6.1 


14.101 


00030956042 


WT 


2009-02-07 


54869.3 


1.7 


16.189 


00090007036 


WT 


2009-02-12 


54874.1 


4.6 


21.056 


00090007037 


WT 


2009-02-22 


54884.5 


4.6 


31.431 


00090007038 


WT 


2009-03-04 


54894.5 


3.9 


41.408 


00090007039 


WT 


2009-03-13 


54904.0 


1.1 


50.912 


00090007040 


WT 


2009-03-24 


54914.8 


4.2 


61.763 


00030956054 


WT 


2009-09-30 


55104.5 


3.3 


251.428 



Table 2 

Burst statistics of magnctars 



Magnetar 


T90 


tr 


*/ r 




pa 


Reference 




(ms) 


(ms) 


(ms) 


(10 14 G) 


00 




SGR 1806-20 


161.8 






21 


7.6 


Goeiis et al. (2001) 


SGR 1900+14 


93.4 






7.3 


8.0 


Gogiis et al. (2001) 


IE 2259+586 


99.31 


2.43 


13.21 1.35 


0.59 


7.0 


Gavriil et al. (2004) 


IE 1547-5408 


305 


39 


66 0.17 


2.2 


2.1 


this work 



"From the McGill AXP/SGR catalogue at http://www.physics.mcgill.ca/~pulsar/magnetar/main.html 
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Hernquist, L. 2004, ApJ, 605, 378 
Woods, P. M., Kouveliotou, C, Gavriil, F. P., Kaspi, M., V., 

Roberts, M. S. E., Ibrahim, A., Markwardt, C. B., Swank, 

J. H., & Finger, M. H. 2005, ApJ, 629, 985 
Woods, P. M. & Thompson, C. 2006, in Compact Stellar X-ray 

Sources, ed. W. H. G. Lewin k, M. van der Klis (UK: 

Cambridge University Press) 
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Figure 1. Properties of the persistent emission of IE 1547—5408 surrounding the 2009 outburst event. The panels on the left show the 
two observations that preceded the event, on a linear time scale. The right-side panels are the post-event observations and are plotted on a 
logarithmic time scale. The dashed line in the top panel shows the power-law decay fit to the unabsorbed flux. The unabsorbed flux in the 
top panel and the pulsed fraction in the bottom panel are for the 1-10 keV range. The first ten points (the first day) after the BAT trigger 
are contaminated by dust scattering and should thus be regarded with caution. Some of the early PC mode observations had background 
regions that were contaminated by the dust scattering; these are marked with an open triangle. 
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Figure 2. 2-10 keV RMS pulsed flux evolution of IE 1547-5408 determined from Swift XRT and RXTE. The black crosses show the 
RXTE pulsed count rates and the red points are the Swift pulsed count rates, arbitrarily scaled to the RXTE values. The dotted vertical 
lines mark the onsets of the 2008 and 2009 outbursts. 
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Figure 3. 1-10 keV RMS pulsed fraction as a function of 1-10 keV unabsorbed flux. Sec also Ng ct al. (2011). 
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Figure 4. Examples of bursts from IE 1547—5408. Time series are binned on a 1/16 s timescale in the 0.5-10 keV energy range. They 
display a wide range of morphology. 




Figure 5. Top left: Distribution of T90 duration of bursts. Top right: Distribution of burst rise times. Bottom left: Distribution of burst 
fall times. Bottom right: Distribution of t r /tf. In all panels, the solid line is the best- fit log- normal function. 
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Figure 6. Left: Distribution of burst fluences. Right: Distribution of burst peak fluxes. The solid line is a linear fit to the filled circles. 
These distributions are based on burst counts from a 1—10 keV energy band. The open circles are not included in the fits because of reduced 
sensitivity in detecting such bursts. The top axes show the fluence and peak flux in cgs units which are derived from a single conversion 
factor between counts and ergs cm -2 . In using such a factor, the bursts are assumed to have the same spectrum which is an approximation 
as each burst has a slightly different spectrum. 
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Figure 7. (a) Folded profile of the times of the burst peaks, (b) Folded profile of photon counts from cycles of the pulsar that contain 
bursts, (c) Folded profile of photon counts from cycles containing symmetric bursts, (d) Folded profile of photon counts from cycles 
containing slow-fall bursts, (e) Folded profile of photon counts from cycles of the pulsar that do not contain bursts but not including the 
first two observations following the BAT trigger, illustrating the quiescent pulse profile. Profiles are all in the 0.5-10 keV energy range. 
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Figure 8. Power-law index as a function of average absorbed flux over T90 from the spectral fits of the 46 most fluent bursts. The black 
points are the individual bursts and the blue open squares represent the weighted averages of the power-law indices for bursts in logarithmic 
fluence bins. 
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Figure 9. 4-10 keV / 2—4 keV flux hardness, H, as a function of 2-10 kcV flux, F for magnetar outbursts. Both are normalised to their 
quiescent values, H q , F q . For IE 1547—5408, the spectral fits used to determine the hardnesses and fluxes were from this work, and the 
first day of the 2009 outbur st has been omitted due to contamination by the dust-scattered emissi on. For XTE J18 10— 197 the spectral 
parameters were tak en fromlGotthelf & Halpernl <|2007T ). For IE 1048.1-5937 they were taken from lTam et~aTl l[200Sft . For IE 2259+586 
they were taken from lZhu et al.l 1 120081) . Uncertainties on the hardness ratios are shown only for IE 1547—5408; determining those for other 
sources is difficult from the literature, however they are likely comparable to the scatter. 



